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Abstract 

The purpose of this paper is to investigate the scattering by a nonlinear crystal 
whose depth is about the wavelength of the impinging field. More precisely, an 
infinite nonlinear slab is illuminated by an incident field which is the sum of three 
plane waves of the same frequency, but with different propagation vectors and 
amplitudes, in such a way that the resulting incident field is periodic. Moreover, 
the height of the slab is of the same order of the wavelength, and therefore the 
so-called slowly varying envelope approximation cannot be used. In our approach 
we take into account some retroactions of the scattered fields between them (for 
instance, we do not use the nondepletion of the pump beam). As a result, a system 
of coupled nonlinear partial differential equations has to be solved. To do this, 
the finite element method (FEM) associated with perfectly matched layers is well 
suited. Nevertheless, when using the FEM, the sources have to be located in the 
meshed area, which is of course impossible when dealing with plane waves. To get 
round this difficulty, the real incident field is simulated by a virtual field emitted by 
an appropriate antenna located in the meshed domain and lying above the obstacle 
(here the slab). 

1 Introduction 

The development of photonic science in nanotechnologies requires an always increas- 
ing control of light. Surface-phenomena, metamaterials or the use of nonlinear optics 
are very efficient ways to do this. In this paper, we combine the last two options: an 
electromagnetic field induces, in a homogeneous nonlinear medium, a periodicity of 
period close to the considered wavelength. A precise description of the system is given 
in the next section. 

The majority of work in nonlinear optics applies to the propagation of a wave in 
a nonlinear medium (see JT] for a review), that is, the wavelength is small compared 
to the length of the path of light in the nonlinear medium. In this case, the paraxial 
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approximation is often used and the equations obtained are of parabolic kinds. For 
example, when studying the propagation of a soliton in an optical fiber (say, oriented 
along the z-axis), one usually restrict the problem in the (z,t)-plane, neglecting the 
transversal effects. This leads to equations of the nonlinear Schrodinger type. Some 
works have been done to determine the limits of this approach (|2][3]|4]). This justi- 
fies the studies done directly with Maxwell's equations, as in J5] [6] , where transversal 
effects are taken into account. 

On the other side, when the wavelength is far larger than the obstacle, a mean-field 
approximation is used ( (Tj [8j ) . The purpose of this paper is to stand between these two 
states, in the realm of resonance, where rough approximations cannot be applied. To 
this end, we use the FEM: the precision in the change of apparent permittivity inside 
the slab (i.e., the inhomogeneity of the sources, resonating with different frequencies, 
of the total field) is only dictated by the memory size of the computer. 

The study of nonlinear optics in periodic material is of course not new (ll9l [lOl . 
see also fiD ). It is known that Kerr photonic crystals allow to produce systems with 
hysteresis (171 [121 ). or with a transmission range that depends on the amplitude of the 
fields ( irT3lfT4l ). Peak of transmission can also appear in the gap due to the nonlinearity 

(Hskhi). 

When considering harmonic generation due to a nonlinear periodic structure, the 
phenomena that appear in photonic crystals are still richer, for the frequency compo- 
nents of the fields can have completely different behaviors ( lfT71[T8l[T9ll20l ). 

Finally, this paper also aims at exposing an application of the companion paper 
IF2TI . that gives a new route to obtain propagation equations in nonlinear optics. 

2 Set up of the problem 

2.1 Description of the system 

The choice of the test structure has been dictated by two guidelines: it must be at the 
same time simple enough to lead to tractable numerical models and possibly feasible 
experiments, and enhance the nonlinear effects both qualitatively and quantitatively. 
For this, we propose the following experiment: let three plane waves impinge on a 
slab, made up of a nonlinear and non-centro-symmetric crystal. One wave is directed 
normally with respect to the slab, the other two are symmetrically oriented with respect 
to it and have the same amplitude (the analytical expression of the incident field is 
given in subsection 13. 11 1. In this way, the problem is periodic, along one direction 
that we call the x-axis. This incident field create an optical lattice ([22, 23 j), as seen 
in the figure Q] Hence, the scattered field presents several orders of reflection and 
of transmission. Furthermore, the scattered field oscillate at several pulsations, due to 
harmonic generation in the crystal. The range of ratios wavelength/period of the crystal 
offers a variety of response for the different harmonics. In particular, the smaller the 
wavelength, the more order excited. This implies that there are more scattered angles 
for the fields oscillating at the generated harmonics than for the one oscillating at the 
incident frequency. Moreover, and we find it spectacular, the energy that flows along 
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each direction of the crystal is not at all a monotonous function of the incident intensity. 
This offers new ways to control the directions along which the higher harmonics escape 
from the nonlinear medium. 

Let Xj be the wavelength of the incident field. The slab we choose had a thickness 
of 7.06 A/. In this range, the slowly varying envelop is not a good approximation, 
so we have to solve the complete set of Maxwell's equations. But, as shown in 11211 . 
Maxwell's equations lead to an infinite set of nonlinear and coupled equations. The 
article just mentioned is mainly devoted to the elaboration of this system, and then its 
truncation, necessary for a numerical study. Some results are briefly reminded in the 
following subsection; the reader is reported to ED for their justification. 

2.2 What kind of equations has to be solved? 

We consider a spatially local, nonbianisotropic, stationary, magnetically linear, smooth 
medium (see [21] for a definition; the idea behind a smooth medium is that we can 
make a Taylor expansion of the electric induction field in function of the electric field). 
For simplicity, the electric susceptibility tensor x are neglected when n > 3. The 
nonlinear effects we want to consider are mainly the second and third harmonic gener- 
ations. Hence the nonharmonic processes (like Raman or Brillioun scatterings) are not 
taken into account. 

We write the electric vector, at a point s and a time t as 

E(s,i) = 2SRe{E J (s)e- l ^ t +E 2 (s)e^ aw/ * +E 3 (s)e- 3l ^ t }. 

Note that this expression implies that we neglect the cascading effect that create har- 
monics higher than the third one; also, as it is often the case, the static component is 
supposed to vanish. 

To write the propagation equations, we need to introduce two notations: the first 
one is the operator M 1 ™, that, when applied to the p-th component E p of the electric 
field, gives the equation satisfied by E p in a linear medium; the second one is a short 
way of writing the interactions between several components of the total electric field. 
For example, [Ei,EiJ := x^ w/)EiEi, describes the interaction of Ei with 
itself - this term is the source of the second harmonic generation -, |_E-i, Ei, EjJ := 
X^ (— Wj, uji, Wf)E_iEiEi, is the term that describes the interaction between E_i, 
Ei and Ei. Since E_ p = E p , this is precisely the optical Kerr effect. More formally, 
we have 

M p m (E p ) := -e^V x (/^(pw/JV x E p ) + {pu^e^ipui)^ 

and 

LE Pl) -- - ,E P J := x (n) 0iwj,-- - ,p n w 7 )E Pl •••E p „ 

where pi G Z. Note that this term contributes to the component of the n-th order 
polarization oscillating at the pulsation (pi + ■ ■ ■ + p n )vi 

Now, as said above, we want to study the second and third harmonic generations. 
Their respective sources are |_Ei,EiJ for E 2 , and |_Ei,E 2 J and |_Ei,Ei,EiJ forE 3 . 
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Figure 1 : The geometry of the system is a slab, invariant along the horizontal direction, 
called the x axis, and along the direction perpendicular to the plane of the sheet of 
paper, called the z axis. Made up of /3 Barium Borate (BBO), it is illuminated by 
three plane waves: one with a normal incidence and an amplitude Ao, and two with 
a symmetric oblique incidence and amplitude A\. The interferences between these 
different incident waves oscillating at cjj, induces a periodic structure, equivalent to 
a diffraction grating or a finite photonic crystal. This is seen in color inside the slab, 
where the pattern of the intensity field ||EjJ| of the incident waves is shown. This 
nonlinear induced grating generate some harmonics. The number of propagating order 
increases with the pulsation. (Three periods are presented in this picture, though the 
simulation runs with one period only.) 
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Note that we consider here a cascading effect: Ei and E2 interact to create E3. In fact, 
we give a set of equations that allows to consider a lossless process ([21, 24|): 

M in (Ei) (la) 
+ w|(2[E_2,E 8 J + 2[E_i,E 2 J 

+ 6[E_ 3 , Ei, E 3 J + 3|_E_i, E_i, E 3 J + 3[E_ 3 , E 2 , E 2 J 
+ 6LE_ 2 ,Ei,E 2 J +3LE_i,Ei,Eij) =s(E\), 

X 2 m (E 2 ) (lb) 
+ (2 W/ ) 2 (2LE_i,E 3 J + [Ei, Ei J 
+ 6[E_ 3 , E 2 , E 3 J + 6[E_ 2 , Ei, E 3 J + 3[E_ 2 , E 2 , E 2 J 
+ 6LE_i,Ei,Eij) =0, 

M m (E 3 ) (1c) 
+ (3w/) 2 (2LEi,E 2 J 

+ 3LE_ 3 ,E 3 ,E 3 J +6LE_ 2) E 2 ,E 3 J +6LE_i,Ei,E 3 J 
+ 3LE_i,E 2 ,E 2 J + LEi,Ei,Eij) =0. 

Note that some other interesting effects are taken into account: the optical Kerr ef- 
fect ( [E 1 , E 1 , E _ 1 J ), the depletion of the pump beam ( [E 2 , E _ 1 J ) . On the other hand, 
some terms are really weak (typically [E 3 , E 3 , E_ 3 J as compared to [_Ei, Ei, EiJ) 
and do not have important consequences on the fields. The presence of these interac- 
tions is only due to the energy conservation property, which is a good test from the 
numerical point of view. Finally, s(E^) is the external source of Ei; since in our case 
the slab is illuminated by plane waves, s(E\) vanishes. 

3 Some Numerical Results 

3.1 From a practical point of view 

To solve the system of equations we posed, numerical methods have to be used. The 
nonlinearity is tackled by a Newton-Raphson scheme and, at each step, the finite el- 
ement method was chosen for its ability to treat the nonhomogeneous sources, in- 
duced by the nonlinearity, of each equation. The incident wave is monochromatic, 
E l (s, t) — 2Ke{Ej (s)e _lW/ *}, and to simplify the comparison with an experimental 
setup, we fix uj such that the associated wavelength is A/ = 1064 nm, which corre- 
sponds to a Nd:YAG laser. 
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Only two-dimensional problems are addressed here. We recall that the physical 
system is invariant along the x and the z axis. The TH polarization is chosen, that is, 
a function u % : M. 2 — > C satisfies, in Cartesian coordinates, 

E\(x,y,z) = u\x,y)z. 

Now a key step is described: the crystal we choose is the (3 Barium Borate (from now 
on denoted as BBO); it is known that there exists one orientation of the crystal that 
guarantees that the total field has the same polarization as the incident field: 

E p (x,y, z) = u p (x,y)z 

for some functions u p : R 2 — > C. We choose this orientation, so that the p-th compo- 
nent of the electric field is described with these scalar functions. 

Finally, the propagation equation system (0 being given in term of the total electric 
field, we use a virtual antenna to simulate the incident field ([25]). Perfectly matched 
layers are used to impose outgoing wave conditions above and below the slatQ, for u\, 
U2 and U3. For the reader's convenience, we give the electric current that simulates the 
incident electric field: 

j(z,y,t)=2Ke{j(a!,y) e"*"'}*, 

with 

j(x, y) = —{ kA e- tky + 2f3A l e" if3y cos(ax)}S(y - h), 

to be evaluated at h, where h is the height of the thread on which the current flows, 
a := k • x = k sin(8) and j3 := k • y = k cos(6*) (k := ||k|| — 2tt/X] is common to the 
three incident waves and 8 is the obliquity angle of the waves with the amplitude Ai), 
and S is the Dirac distribution. The expression of the incident field is then 

u l (x, y) = A e~ tky + 2A ie ' ll3y cos(aa;). 

3.2 Scattering Of Three Waves On A Slab 

We come now to the results of the simulation of the experiment described in the figure 
Q] This system exhibits a non trivial and quite complex behavior because of the induced 
diffraction grating. 

The numerical value of the relevant susceptibility components are given in the table 
[T]We note that they do not depend on the frequency. 

We are interested in the directions where the fields scatter. Out of the slab, for each 
frequency, the electric field satisfies a Helmholtz equation. We write its propagating 
solution as 

'From the physical point of view, only the scattered part of u± satisfies an outgoing wave condition. But 
the principal of the virtual antenna is precisely to find a source, located in the meshed domain, that generates, 
in the neighborhood of the obstacle, exactly the incident field. All the details are given in f 25l . 
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Table 1 : The wavelength, in nm, the refraction index of the extraordinary axis, and the 
values of the relevant (in that article) x' 1 ), an d X > respectively without unit, in 
m.V^ 1 and in m 2 .V~ 2 . 



+ 6 W e <^(^-(^) 2 ) 1 N}£, 

(r) (t) 

where b p ^ (resp. b p ; n ) denotes the coefficient of the n-th order of the reflected (resp. 
transmitted) wave at pulsation pujj, and U p :— {n G Z : k 2 — (^p) 2 > 0}. Since 
k p = puji/c (where c is the light velocity), the higher the harmonic p, the higher 
the modulus of the wave vector k p ; in other words, the size of U p increases with p. 
This means that there are more scattered angles for higher harmonics than for the field 
oscillating at the fundamental frequency. In our particular case, whereas only the orders 
— 1, and 1 are present in Ei, the third harmonic contains orders from —3 to 3. These 
ranges depend on the wavelength and the inclination 9 of the incident beams; had we 
chosen a larger d, the propagating orders allowed for the first harmonic would have 
ranged from, say, —2 to 2. 

The second thing we note is that the amplitude of the n-th order of the p-th har- 

(r) (t) 

monic (that is b Pt n or b p ' n ) is not monotonic in the amplitude of the incident field, as is 
seen in the figures |2]|4] (due to the reflection symmetry along the y-axis of the system, 
and b p r l n = b [ ;} n ). In the fi gures are represented the efficiency of the 
transmitted waves - the reflected ones show similar behaviors - defined by: 

\h (a) \ 

V P,n ■ „ , (a) 

l^nGUp \ u P, n \ 

for a G {r, t). It thus gives the part of the intensity of the reflected or transmitted wave 
at pulsation pui that escapes along the n-th order. 



4 Concluding Remarks 

We present numerical evidences that the scattering of several waves on a nonlinear 
slab, whose thickness is of the order of the wavelength, shows rich phenomena. This 
system is studied through a set of equations obtained from a rigorous method, given in 
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(a) (b) e« 

Figure 2: The efficiency of the transmitted propagating waves at the frequency w/. As 
is easily understood, when A\ tends to zero, all the light escape perpendicularly to the 
slab, i.e., tends to one and e± \ vanishes. 




Figure 3: The efficiency of the transmitted propagating waves at the frequency 2wj. 
Schematically, when A Q > A x (resp. A ~ A lt A Q < A x ), the th (resp. I s *, 2 nd ) 
order is favored. That \ is peaked around Aq = Ai was a surprise for us. 
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a companion paper. The response of the induced grating is far from being monotonous 
in the amplitudes of the pump wave. We now look for an experimental test. 
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